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FOREWORD 


The  design  and  development  of  liquid  filled  acoustic  lenses  necessitated  a 
theoretical  method  which  can  accurately  predict  the  transmission  coefficient  of 
an  ac  istic  wave  propagating  through  layered  media  at  various  angles  to  the 
normal.  This  report  describes  such  a  method  for  determining  the  transmission 
coefficients  and  demonstrates  its  accuracy  against  available  published  data. 
Funding  for  this  effort  was  provided  by  the  Naval  Sea  Systems  Command,  Task  No. 
S0266001/U13CA. 

The  author  would  like  to  thank  Dr.  Bruce  Hartmann  for  his  guidance  in  this 
endeavor,  Mr.  Paul  Huber  for  his  assistance  in  verifying  the  accuracy  of  the 
method  against  published  data,  and  Mr.  V.  Rust,  Head,  Electrical  Design  Branch, 
for  reviewing  this  report. 
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INTRODUCTION 

Many  materials  are  used  in  underwater  acoustic  systems  for  items  such  as 
transducer  windows,  acoustic  concentrators,  acoustic  lens,  acoustic  tank 
linings,  and  reflectors.  Their  application  depends  strongly  on  their 
characteristic  properties  in  transmission,  reflection,  and  absorption.  The 
continuous  development  in  sonar  domes  and  acoustic  lens  necessitated  the 
investigation  of  ultra-sonic  wave  through  single  and  multilayered  thin  plastic 
materials.  Until  recent  years  there  has  been  limited  experimental  data  and 
theoretical  results.  Often  the  transmission  and  reflection  coefficients  are 
measured  instead  of  calculated.  To  date,  there  seems  to  be  several  methods 
available  for  theoretically  approximating  the  transmission  characteristics  of 
the  acoustic  materials.  This  report  presents  one  of  the  methods  and  its 
calculated  results. 

One  of  the  many  applications  of  acoustic  materials  in  underwater  acoustics 
is  in  the  construction  of  acoustic  lens.  In  particular,  this  program  has  an 
interest  in  a  spherical  liquid  filled  lens  as  diagrammed  in  Figure  1.  A 
liquid-filled  spherical  lens  basically  consists  of  a  thin  plastic  shell  of  the 
desired  transmittion  coefficient  and  a  liquid  of  the  desired  index  of  refraction 
(h2^ •  It  has  the  properties  of  converging  the  acoustic  rays  toward  a  focal 
point  along  the  lens  axis  as  shown  in  Figure  1.  Its  focal  distance,  f,  from 
center  of  the  spherical  lens  is  directly  related  to  the  index  of  refraction  of 
the  lens  liquid. 

To  select  the  proper  acoustic  window  for  an  acoustic  lens,  it  is  necessary 
to  select  the  transmission  characteristics  of  alternative  materials.  These 
characteristics  are  evaluated  on  the  basis  of  determining  the  percentage  of 
sound  transmitted  through  and  reflected  from  the  sample  material  immersed  in  a 
liquid  media.  Through  the  transmission  characteristic  of  a  material,  one  can 
predict  the  approximate  aperture  of  an  acoustic  lens.  For  a  desired  aperture, 
material  of  certain  transmission  characteristic  will  be  selected.  Having 
selected  a  desired  aperture  one  would  have  defined  the  desired  operational 
characteristic  such  as  beam  width  of  the  acoustic  lens. 

The  problem  of  transmission  and  reflection  of  ultrasonic  waves  through 
layered  media  has  been  treated  by  many  authors.  The  theoretical  treatment  dates 
back  to  W.  T.  Thomson!  in  1950.  In  nis  paper,  he  treated  the  transmission  and 
reflection  problem  as  a  boundary  value  problem.  At  the  boundary,  the  stress  and 


!Thomson.  W.,  "Transmission  of  Elastic  Waves  Through  a  Stratified  Medium, 
J.  Appl.  Phys.  21,  39-93,  1950. 
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particle  velocity  are  assumed  to  be  continuous.  Brekhovskikh’s^  treatment 
paralleled  that  of  Thompson's  and  yields  a  similar  expression  which  is  valid  only 
to  limited  cases:  i.e.,  single-layer  problem  as  indicated  by  Mr.  D.  Folds.3 
In  1953,  Haskell^  correctly  formulated  an  expression  which  is  applicable  to 
multilayered  cases.  Since  our  interest  is  a  single  solid  layer  immersed  in  a 
liquid  media  and  because  Brekhovskikh's  solution  is  a  widely  used  reference  by 
many  researchers,  we  have  employed  his  expression  for  this  analysis. 


2Brekhovskikh,  L.  M.,  "Wave  in  Layered  Media."  Academic,  New  York,  I960. 


3 Folds,  D.  L.  and  Loggins 
Waves  in  Layered  Media," 


C.  D.,  "Transmission  and  Reflection  of  Ultrasonic 
J.  Appl.  Soc.  Am.  62,  1102-1109,  1977- 


^Haskell,  N.,  "The  Dispersion  of  Surface  Waves  in  Multilayered  Media,"  Bull. 
Seismol.  Soc.  Am.  43,  17-34,  1953* 


9/10 


THEORY 


Paralleling  Brekhovskikh ’a  development  with  Mr.  P.  Huber's  assistance  the 
same  expressions  were  derived  for  a  single  solid  layer  surrounded  by  a  liquid 
media  as  D.  Folds  and  Barnard,  et  al5  had  indicated  in  their  papers.  The 
transmission  (t)  and  reflection  (r)  coefficient  expressions  for  a  single  solid 
layer  surrounded  by  liquid  media  as  diagrammed  in  Figure  2  are  given  by: 


2  N  Z 


2  2 

M  (Zx  +  Z3)  +  i(N  -  M  )  Z1  +  Z3 


M(Z  -  Z  )  +  i(N2  -  M2)  Z  -  Z 

R  =  --A-.-3 _ 1  3 

2  2 

M  (Zx  +  Z3)  +  i(N  -  M  )  Z1  +  Z3 

where  M,  N,  and  Z's  are  expressed  a3: 

Z  Z 

M  =  — L.  Cos2  2y  Cot  P  +  — il-  •  Sin2  2  y  •  Cot  Q 
Z1  J2  Z1  2 


Z  Cos2  2y  Z  Sin2  2y 

N  =  _2 _ 2  «.  2t  . _ 2 

Z±  Sin  P  Z^  Sin  Q 


and 


P  o 

Z  =  i-  1 
1  Cos 


^Barnard,  G. ,  Bardin,  J.  L.  and  Whiteley,  J.  W. ,  "Acoustic  Reflection  and 
Transmission  Characteristics  for  Thin  Plates."  J.  Acoustic  Soc.  Am.  57- 

577-584,  1975. 
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r'  'T*''  1  *4*  am 


Z 

2 


P  c 


2  2 


Cos  0 


2 


P  c 

Z  =  -  2 . 3._ 

3  Cos  9  j 


2t  Cos  7 2 

where  the  subscripts  1,  2,  and  3  indicate  the  medium  of  interest  while  2t 
represents  the  shear  (transverse)  mode  of  propagation  in  solid  layer,  media  2. 
Also,  p's  are  the  densities  and  c's  and  b's  are  the  velocities  of  the  respective 
media;  P  and  Q  are  two  "dummy"  variables  which  are  related  to  the  propagation 
constant  for  the  longitudinal  wave,  *  ,  shear  wave,  0,  and  thickness,  d,  of  the 
sample  material 


P  =  ad 
Q  =  /3d 

Since  a  represents  the  z-component  of  the  propagational  constant  k2  of  the 
longitudinal  wave  while  0  is  the  like  quantity  K£  of  the  shear  wave,  thus  a 
and  0  can  be  expressed  as 

a  =  -k2  Cos  0  2 

0  =  -K2  Cos  >2 

where  #2  and  ^2  are  related  to  the  incident  angle  #3  by  the  Snell's  Law 

k3  sin  9 3  =  k2  sin  9 2  =  k]_  sin  0^  =  K2  3in  y2 

9 3,  62 1  9 1,  and  y2  are  measured  with  respect  to  the  normal  to  the  plane 

of  the  solid  layer  and  these  quantities  may  be  complex.  Thus,  the  quantities  P 

an  Q  can  be  rewritten  as; 

P  =  -k2d  Cos  (?2 
Q  =  -K2d  Cos  y 2 

Note  that  k2  and  K2  are  wave  numbers  in  the  longitudinal  and  shear  mode 
which  are  given  respectively  by 

k  =  Ji 
2  c2 
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W  =  2r  •  f 

where  03  and  b2  are  respectively  the  longitudinal  and  shear  sound  velocities 
in  the  solid  layer  of  thickness,  d,  and  f  is  the  frequency  of  interest. 

Since  absorption  is  present  in  most  solid  materials  and  it  is  related  to 
the  sound  velocity  c  and  b  in  the  solid  by  the  expression  below 


c 

c  =  _ 2 _  EXP  (i  tan  -1  <*>  ) 

2  1/2  c 

(1  ♦  *c  ) 


where 


a  and  0  are  redefined  as  the  absorption  constants  expressed  in  neper/ft.  The 
absorption  can  be  related  to  the  attenuation  by 

a  =  a  /  8.6858 

where  a  is  the  attenuation  in  db/ft.  Note  that  0£  and  bg  are  now  complex 
quantities.  It  is  indeed  difficult  to  image  physically  that  the  3ound  velocity 
and  reflected  angles  in  the  solid  layer  as  complex  quantities. 
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EXPERIMENTAL  DATA 

The  experimental  test  set  up  at  NAVSURFWPNCEN  acoustical  facility  is 
basically  that  of  References  3#  5,  and  6,  except  for  the  collimator  lens  which 
was  used  as  diagrammed  in  Figure  3.  The  collimator  was  used  to  insure  that  the 
impinging  wave  at  the  sample  plate  would  be  a  plane  wave  (far-field  effect). 

The  receive  hydrophone  E-27  is  located  close  to  the  sample  plate  to  avoid  edge 
effects;  yet  far  enough  to  avoid  standing  waves  between  the  hydrophone  and  the 
sample  plate.  The  test  sample  plate  was  a  24-inch  square  and  hydrophone 
locations  are  as  shown  in  Figure  4. 

Since  only  limited  experimental  measurements  were  made  at  NSWC  and  the 
sound  velocities  and  test  material  are  not  known,  experimental  data  used  to 
verify  calculated  results  were  basically  from  References  3  and  6.  Figure  5 
to  Figure  8  show  results  from  Reference  3  versus  NSWC  calculated  results. 

Figure  9  to  Figure  16  show  the  data  from  Reference  3  versus  NSWC  calculated 
results.  The  material  parameters  from  References  3  and  7  used  in  the  theoretical 
calculations  are  listed  in  Table  1. 


3see  footnote  3  on  page  9 . 
5see  footnote  5  on  page  11  . 


^Nichols,  R.  H. ,  "Acoustic  Properties  of  Rho-C  Rubber  and  ABS  in  the  Frequency 
Range  100  kHz  -  2  MHz,"  J.  Acoustic  Soc.  Am.  54,  1763-1765,  1973. 


^Hartmann,  B.  and  Jarznyski,  J. .  "Immersion  Apparatus  for  Ultrasonic 
Measurements  in  Ployraers,"  J.  Acoustic  Soc.  .Am.  56,  1479-1477,  19711, 


15 


NSWC  TR  81-197 


RESULTS 


Brekhovskikh ' s  expression  was  used  in  the  theoretical  calculation  to 
predict  the  transmission  coefficients.  To  better  estimate  the  coefficients,  two 
assumptions  were  made.  First,  the  velocities  c  and  b  were  allowed  to  be 
complex;  thus,  the  wave  numbers  in  the  solid  layer  k£  and  K2  are  also 
complex.  These  complex  quantities  were  used  throughout  the  calculations. 
Secondly,  instead  of  assuming  the  attenuations  in  the  longitudinal  and  shear 
modes,  a  and  0  respectively  are  the  same,  they  are  assumed  to  have  the  relation 
of 


_a_  -  20  / b^ 

\  xs  \<y 

Since  X^  and  X3,  the  respective  wavelength,  are  defined  as 


X 


s 


thus,  the  above  can  be  rewritten  as 

*  •  23 (£) 


and  0  can  be  expressed  as 


0 


Throughout  the  development  of  the  method  of  approximation  it  has  been 
observed  that  the  correct  sound  velocity  measurements  are  very  important.  For 
example,  in  determining  the  theoretical  transmission  coefficient  for  Absonic-A 
at  19.1°C  of  Reference  5,  there  was  no  difficulty  arriving  at  the  same  data  as 
the  measured  data  when  measured  parameters  of  Reference  5  were  used. 

These  results  were  shown  in  Figures  5  through  8.  However,  if  other  published 


5see  footnote  5  on  page  11. 
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values  were  used  for  the  Absonlc -A  sound  velocities,  the  calculated  results  were 
very  different  from  the  measured  data.  In  the  case  of  Reference  3.  using  the 
values  in  Table  1,  the  calculated  transmission  coefficient  was  a  good 
approximation  of  the  published  measured  data  for  a  single  solid  layer.  This  is 
shown  in  Figure  9  for  LDPE;  Figures  10  through  13  for  plexiglas;  and  Figures  14 
and  15  for  ABS. 

It  is  evident  from  this  theoretical  investigation  that  the  material 
transmission  coefficient  is  dependent  on  frequency,  temperature  and  thickness  of 
solid  layer.  For  example,  Figures  10  through  13  show  the  effect  on  the 
transmission  coefficient  of  a  0.1  inch  thick  plexiglas  plate  as  a  function  of 
frequency.  A  similar  effect  is  shown  on  Figures  5  through  8  and  Figures  10 
through  18  for  Absonic-A  (ABS).  This  frequency  effect  is  observed  both 
theoretically  and  experimentally  as  shown  in  these  plots.  Since  sound  velocity 
changes  as  a  function  of  temperature,  one  could  expect  some  change  in  the 
transmission  coefficient.  In  the  case  of  the  materials  considered,  the 
temperature  effect  is  amplified  greatly  because  sound  velocity  increases  with 
increase  of  temperature  in  the  liquid  (water)  while  the  contrary  exists  in  the 
solid  layer.  The  observed  results  are  shown  in  Figures  19  though  23  for 
low-density  polyethylene  (LDPE).  However,  material  with  low  thermal 
coefficients  such  as  ABS  has  very  little  effect  due  to  temperature  change  as 
shown  in  Figures  24  through  26.  Thus  careful  selection  of  material  could  reduce 
the  temperature  effect  on  the  transmission  coefficient.  Lastly,  thickness  of 
the  shell  also  has  sane  effect  on  the  transmission  coefficient  as  shown  in 
Figures  27  through  36.  From  these  plots  one  could  observe  the  changes  on 
transmission  coefficient  as  a  functon  of  thickness. 

Aperture  shading  is  well  known  and  has  been  used  successfully  in  optics  and 
underwater  acoustics.  In  general,  aperature  shading  is  used  to  improve  the 
directivity  response  of  an  acoustic  system.  The  transmission  coefficient  on  an 
acoustic  lens  shell  has  the  same  effect  as  that  of  a  shading  function.  Thus,  a 
lens  of  a  certain  operating  characteristic  would  require  a  certain  desired 
transmission  coefficient  on  its  shell.  Througn  careful  consideration  of  the 
material,  thickness,  frequency,  and  operating  temperature  range,  a  lens  shell 
can  be  selected  to  exhibit  a  desired  transmission  coefficient  with  aperture. 


3see  footnote  3  on  page  9. 
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CONCLUSION 

It  is  evident  from  this  investigation  that  the  acoustic  transmission  and 
reflection  coefficient  of  a  thin  plate  can  be  accurately  predicted  by  employing 
the  Brekhovskinkh' s  expressions  in  conjunction  vith  an  accurately  determined  sound 
velocity  in  the  thin  plate  material.  Thus,  this  technique  can  be  useful  in 
assisting  the  selection  of  acoustic  material  for  acoustic  lens  of  the  operating 
characteristic  desired  for  the  acoustic  system. 


I  FIGURE  S  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  100  kHz  FOR  ABS  MATERIAL 


FIGURE  6  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  300  KHz  FOR  ABS  MATERIAL 
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FIGURE  7  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  500  kHz  FOR  ABS  MATERIAL 


FIGURE  8  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  1  mHz  FOR  ABS  MATERIAL 


FIGURE  9  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE 

AT  500  kH*  FOR  LOW  DENSITY  POLYETHYLENE  MATERIAL 


FIGURE  10  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  100  kHz  FOR  PLEXIGLAS  MATERIAL 


FIGURE  11  TRANSMISSION  LOSS  AS  FACTION  OF  INCIDENT  ANGLE  AT  200  kHz  FOR  PLEXIGLAS  MATERIAL 


FIGURE  15  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  200  kHz  FOR  ABS  MATERIAL 


FIGURE  16  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  300  kHi  FOR  ABS  MATERIAL 
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FIGURE  17  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  400  kHz  FOR  ABS  MATERIAL 
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FIGURE  18  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  500  kHz  FOR  ABS  MATERIAL 


FIGURE  19  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE 
AT  20° C  FOR  LOW  DENSITY  POLYETHYLENE  MATERIAL 
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FIGURE  22  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE 
AT  5°C  for  low-density  polyethylene  material 
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FIGURE  23  TRANSMISSION  LOSS  AS  FUNCT  ION  OF  INCIDENT  ANGLE 
AT  0°C  FOR  LOW-DENSITY  POLYETHYLENE  MATERIAL 


FIGURE  24  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  20°C  FOR  ABS  MATERIAL 


FIGURE  25  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  10. 5° C  FOR  ABS  MATERIAL 


FIGURE  26  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  AT  0°C  FOR  ABS  MATERIAL 


FIGURE  27  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  FOR  .05  INCH-THICK  ABS  MATERIAL 


FIGURE  30  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  FOR  .2  INCH-THICK  ABS  MATERIAL 


FIGURE  31  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  FOR. 25  INCH-THICK  ABS  MATERIAL 
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FIGURE  32  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  FOR  .3  INCH-THICK  ABS  MATERIAL 
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FIGURE  33  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  FOR  .4  INCH-THICK  ABS  MATERIAL 
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FIGURE  34  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  FOR  .5  INCH  THICK  ABS  MATERIAL 


FIGURE  35  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  FOR  .75  INCH-THICK  ABS  MATERIAL 
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FIGURE  36  TRANSMISSION  LOSS  AS  FUNCTION  OF  INCIDENT  ANGLE  FOR  1  INCH- THICK  ABS  MATERIAL 


NSWC  TR  81-197 


so 

0) 

,o|-p 

•  n 


M 

« 

ol-p 
n3|TS  o 

•  « 


>> 

CM 

$  CO 

on 

cm 

vO 

ON 

non  o 

fH 

O' 

as 

«H 

cn  g  • 

GO  -H 

• 

fH 

• 

o 

• 

o 

CM 

c 

o  tj 

•r-l  35 
Vi  -p  5* 
d  PiS 
di  3< 

3=  °  ^ 
cfl  cn  dq 

A  -a 
a) 


3  p  o 

p  p<  *— 1 
•H  S- 
Id  O  B 
CON 
Q  P  31 
kJ  cd  -O 


>i  O 

U  <u 
a  •—  « 
«  o  ^ 

323 

4> 

> 


T3  >»  O 
3  -p  ail 

*>  H  Ml 
•H  O  S 
00  O  SI 
C  f"4  * 

3  5 


O  o 

UO  LTi 


o  in 

ao 

CM 


O 

o 

O' 

o 

o 

o 

o 

o 

CM 

o 

on 

m 

^r 

cm 

CM 

00 

rH 

0) 

>»  O 

>»  c 

H  G 

H  O 

•h  a> 

/H  < 

03 

■H  <H 

CO  H 

flj  1 

03 

01  >i 

C  >N 

•H  O 

fH 

c  £ 

0)  .C 

im  *H 

ao 

01 

■°  ^ 

c 

o  c 

•H 

03  0 

1 

0 

H  O 

fli  0) 

X 

0) 

1  >* 

36  h 

rH 

4-» 

X  .O 

rH 

0  o 

•H  0 

<U 

£  <c 

a- 

-4  a 

x  a 

H 

footnote  7  on  page  15 
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